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Abstract 

Coherent beam-beam tune shift of unsymmetrical beam-beam interactions was studied experi- 
mentally and numerically in HERA where the lepton beam has a very large beam-beam parameter 
(up to S^y = 0.272). Unlike the symmetrical case of beam-beam interactions, the ratio of the coher- 
ent and incoherent beam-beam tune shift in this unsymmetrical case of beam-beam interactions 
was found to decrease monotonically with increase of the beam-beam parameter. The results of 
self-consistent beam-beam simulation, the linearized Vlasov equation, and the rigid-beam model 
were compared with the experimental measurement. It was found that the coherent beam-beam 
tune shifts measured in the experiment and calculated in the simulation agree remarkably well but 
they are much smaller than those calculated by the linearized Vlasov equation with the single-mode 
approximation or the rigid-beam model. The study indicated that the single-mode approximation 
in the linearization of Vlasov equation is not valid in the case of unsymmetrical beam-beam in- 
teractions. The rigid-beam model is valid only with a small beam-beam parameter in the case of 
unsymmetrical beam-beam interactions. 
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I. INTRODUCTION 

To achieve a substantial increase of luminosity in a storage-ring collider, limited options include 
increase of bunch currents, reduction of beam sizes at interaction points (IPs), and increase of the 
number of colliding bunches. The first two measures unavoidably increase head-on beam-beam 
forces which could lead to collective (coherent) beam-beam instabilities ^H]. Understanding of 
coherent beam-beam effects especially in the nonlinear regime is therefore of primary importance 
for achieving high luminosity in a storage-ring collider with high-intensity beams. 

To study the coherent beam-beam effect, one important quantity that can be measured exper- 
imentally is the coherent beam-beam tune shift. Without beam-beam interactions and without 
considering nonlinearities in the lattice, the two counter-rotating beams oscillate transversely with 
frequencies that correspond to lattice tunes (betatron tunes without collision) if they deviate from 
close orbits. With beam-beam interactions, the particle distributions of the beams are perturbed 
and evolve with time according to the Vlasov equation T. The dynamics of the beams could 
therefore be complicated by multi-mode oscillations of the beam distributions. If considering only 
the stable oscillation of beam centroids (coherent dipole oscillation), the frequency spectrum of the 
beam-centroid oscillation has two primary frequencies for each degree of freedom of the transverse 
motion. These primary frequencies correspond to the tunes measured during collision. The dif- 
ferences between these measured collision tunes and the lattice tunes are the coherent beam-beam 
tune shifts j^Oll^. Over decades, many studies have been conducted on the relationship between 
the coherent beam-beam tune shift and the beam-beam parameter that measures the strength of 
the beam-beam interaction [3-13]. Two theoretical models, the linearized Vlasov equation fTl and 
the rigid-beam model [^, have been studied extensively for cases of weak beam-beam perturbation 
in which the beam-beam parameter is relatively small. When the two beams have the same or very 
close lattice tunes, the calculation of the coherent beam-beam tune shift based on the linearized 
Vlasov equation with single-mode approximation agrees with beam measurements and computer 
simulations |31|^iniIZ]. The rigid-beam model is inconsistent with the linearized Vlasov equation 
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and was therefore proven to be wrong by beam measurements in this case [HllZl- When the two 
beams have very different lattice tunes, on the other hand, the calculation based on the rigid-beam 
model provides a good agreement with beam measurements [S1|H1^^. In both of these models, the 
equilibrium beam distributions were assumed to be Gaussian distributions for easing the calcula- 
tions. In the case of weak beam-beam perturbation, this assumption is fairly good as the beams 
were observed to stay close to a Gaussian. 

The situation of strong beam-beam perturbations with relatively large beam-beam parameter 
is much more complicated and less understood. When the beam-beam parameter exceeds a thresh- 
old, the beam-beam interaction could induce a chaotic coherent beam-beam instability. After the 
onset of the instability, the closed orbits could become unstable for the beam centroids and two 
beams could develop a spontaneous unstable coherent oscillation When the beam-beam 

parameter is below the beam-beam threshold, the coherent beam oscillation is stable. It is, how- 
ever, not clear whether the linearized Vlasov equation or the rigid-beam model are still valid in the 
regime of strong beam-beam perturbation. As many efforts are being made to further increase the 
beam-beam parameter in upgrades of current and developments of future storage-ring colliders, an 
understanding of the coherent beam-beam tune shift in this regime not only is necessary for the 
interpretation of the tune measurement during operation of colliders with high-intensity beams but 
also could shed light on the onset of the chaotic coherent beam-beam stability. 

To explore the beam-beam effect with a large beam-beam parameter a beam experiment, the 
HERA 2000 beam study, was performed on HERA [Hadron Electron Ring Accelerator at DESY 
(Deutsches Elektronen-Synchrotron) , Hamburg, Germany] in which a 920 GeV proton (p) beam 
and a 27.5 GeV positron (e^) beam collided at two IPs, HI and ZEUS ^H]- The beam-beam 
interaction in HERA is a typical case of unsymmetrical beam-beam interaction as the two beams 
have very different lattice tunes and beam-beam parameters (strong p beam and weak e"*" beam). 
In the experiment, the vertical beam-beam parameter of the e"^ beam was varied from 0.068 to 
0.272 by changing the vertical beta-function of the e"*" beam at two IPs. The emittance of the e"*" 
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beam and the luminosity were measured as functions of the beam-beam parameter. One important 
phenomenon observed in this experiment is that the measured coherent beam-beam tune shifts of 
the e"'" beam are much smaller than those calculated from the rigid-beam model. This is the first 
experimental evidence indicating that the traditional models of the coherent beam-beam tune shift 
are no longer valid in the situation of strong beam-beam perturbations. 

To have a better understanding of the experimental data, we reconstructed the HERA beam 
experiment with a self-consistent beam-beam simulation. Remarkable agreement between the ex- 
periment and the simulation was observed on emittance growth and luminosity reduction. More 
significantly, the computer simulation confirmed the experimental result of very small coherent 
beam-beam tune shifts in this case of a very large beam-beam parameter. To examine the validity 
of the theoretical models for the coherent beam-beam tune shift, the linearized Vlasov equation and 
the rigid-beam model were solved for the HERA experiment. Since the distribution of the e+ beam 
significantly deviated from the Gaussian due to the strong beam-beam interaction, the solutions of 
the linearized Vlasov equation and the rigid-beam model were calculated with the beam distribu- 
tions obtained from the beam-beam simulation instead of assuming Gaussian distributions. It was 
found that for the unsymmetrical beam-beam interaction with a large beam-beam parameter, the 
result of the rigid-beam model is inconsistent with the beam experiment and beam simulation even 
though a more accurate beam distribution was used in the calculation. The coherent beam-beam 
tune shifts calculated from the linearized Vlasov equation with the single-mode approximation were 
also found to be significantly different from the result of the beam experiment and beam simulation 
no mater whether the beam-beam parameter is large or small. As the linearization of the Vlasov 
equation is expected to be valid for at least a small beam-beam parameter, this discrepancy sug- 
gests that the single-mode approximation used for solving the linearized Vlasov equation may not 
be valid in the case of unsymmetrical beam-beam interactions. Unfortunately, without the single- 
mode approximation the linearized Vlasov equation for the problem of beam-beam interactions is 
currently unsolvable computationally due to the unsolved degeneracy problem of a matrix with 
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mode coupling. 

This paper is organized as follows. Section II summarizes the results of the HERA 2000 beam 
experiment. In Section III, the self-consistent beam simulation for the HERA beam experiment is 
discussed. In Section IV, the coherent beam-beam tune shifts calculated by using the rigid-beam 
model or the linearized Vlasov equation are compared with the experiment /simulation results. The 
details of the coherent tune calculation with the theoretical models are presented in Appendix A 
and B. The characteristics of the coherent beam-beam tune shift in the unsymmetrical case of 
beam-beam interactions are discussed in Section V. Section VI contains a summary remark. 

II. HERA 2000 BEAM EXPERIMENT 

In the luminosity upgrade of HERA, the beam-beam parameters of the electron beam have been 
nearly doubled. To examine any possible luminosity reduction due to beam-beam effects, a series 
of beam experiments were performed in HERA jl4l I16j . In the HERA 2000 beam experiment, the 
e"*" beam was used to collide with the p beam at the two IPs and the effect of a large beam-beam 
parameter of the lepton beam was explored by increasing the vertical beta-function {(3e.y) of the 
e"*" beam at the IPs. The vertical beam-beam parameter {E,e,y) of the e"*" beam is related to j3e,y by 

m 

^ 1^ ^ (1) 

27r7e (Tp^y{(Tp^x + 0-p,y) 

where and 7e are the classic radius and Lorentz factor of positron, respectively, and Np is the 
number of protons per bunch. The horizontal and vertical size of the p beam at the IPs are given by 
ap^x and (Tp^y In the experiment, the p beam current (/p) was fixed. Since the beam-beam parameter 
of the p beam is very small, there was little change in the p beam size as it was observed during the 
experiment. The vertical beam-beam parameter ^e,y is therefore linearly proportional to (3e,y in this 
case. During the experiment, after the proton current was filled, ^e.y was increased from 0.068 to 
0.272 as Pe,y was changed from 1.0 to 4.0 m while other lattice parameters were kept as constants. 
It should be noted that with two IPs in the HERA experiment a beam-beam parameter of 0.272 
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is among the highest ever achieved in storage-ring colhders. Table 1 hsts the e"*" beam current 
(/e) and the beam-beam parameters of the e"*" and p beam at Pe,y where the measurement was 
performed. For the p beam, Ip = 90 mA, the beta-function at the IPs were {Pp^x,Pp,y) = (7.0, 0.5) 
m; the lattice tunes were (I'xji'y) = (31.291, 32.297); and the emittance without collision was 
i^Oxj^Oy) = (3.82, 3.18) nm. For the e"*" beam, f3e,x = 2.5 m, {vx^i^y) = (52.169, 52.246), and 
(eoxjEoy) = (32.0, 1.28) nm. Other beam parameters can be found in jl6j . The and p beam 
sizes were not matched during the collision. In the experiment, the emittance of the beam and 
the luminosity were measured as functions of Pe,y at both the IPs. In Figs. 1 and 2, the measured 
emittance and the specific luminosity were plotted, with discrete points, as functions of ^e,j/- For 
each ^e,y where the measurement was performed, two data points correspond to the measurement 
at the two IPs, respectively. The specific luminosity is defined as Cg = NcoiC/{IeIp) where N^oi 
and C are the number of colliding bunches and the luminosity, respectively |14| . As shown in 
Fig. 1, the vertical emittance growth of the e+ beam increases monotonically and smoothly with 
the increase of ^e,y This is the characteristics of the incoherent beam-beam effect in contrast to 
the coherent beam-beam effect of which the emittance growth as a function of the beam-beam 
parameter could experience certain jumps (phase transitions) due to the onset of the coherent 
beam-beam instability (2j. To confirm that the luminosity reduction in Fig. 2 is indeed due to the 
emittance blowup, the luminosity calculated with the measured emittance by using the standard 
formula jl6j is also plotted in the Fig. 2. The agreement between the measured luminosity and the 
calculated luminosity in Fig. 2 shows a consistency in the emittance and luminosity measurement. 

In the experiment, the collision tunes of the e"*" beam were also measured at I3e,y = 4.0 m as Vx = 
52.160 and Uy = 52.233. The coherent beam-beam tune shift of the e"*" beam is therefore only Az^^^ = 
0.009 and Af^ = 0.013 while from the rigid-beam model Ai/^. = 0.016 and Az^^ = 0.042 if both the 
beams are Gaussian ^1]. The measured coherent beam-beam tune shifts in this case are inconsistent 
with the traditional understanding of beam coherent oscillation. Moreover, in the symmetrical case 
of beam-beam interactions, the ratio of the coherent and incoherent beam-beam tune shift has a 
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value approximately ranging from 1.2 for round beams to 1.3 for flat (one-dimensional) beams j^. 
Since the beam-beam parameters of the p beam is very small, the distribution of the p beam in 
the HERA experiment is still very close to a Gaussian. The incoherent beam-beam tune shifts 
that are the tune shifts of a single particle with an infinitesimal betatron oscillation are 2^e,a: = 
0.082 and 2,^e,j/ = 0.544 for the e"^ beam in the horizontal and vertical plane, respectively. The 
ratios of coherent and incoherent beam-beam tune shifts of the beam are Ai'e^x/{'^^e,x) = 0.11 
and Ai/e,j//(2^e,?;) = 0.024 at fie,y = 4.0 m. The coherent beam-beam tune shifts in this case of 
unsymmetrical beam-beam interactions with a large beam-beam parameter are therefore extremely 
small as compared with the incoherent beam-beam tune shifts or compared with the symmetrical 
case of beam-beam interactions. 

III. RECONSTRUCTION OF HERA BEAM EXPERIMENT WITH NUMERICAL 

SIMULATION 

To have a better understanding of the measured data in the HERA experiment, we have recon- 
structed the experiment with a self-consistent beam-beam simulation. This numerical study also 
served as a detailed benchmark of our beam-beam simulation code with the experimental measure- 
ment. In the simulation, the linear HERA lattice with the two IPs was used. The code used in 
this study is an expanded version of [2] that is currently capable of studying beam-beam effects of 
proton or lepton beams with any aspect ratio (ratio between vertical and horizontal beam size). 
The two colliding beams were represented by a million macro-particles with given initial Gaussian 
distributions in transverse phase space. Without beam-beam interactions, the initial beam distri- 
bution used in the simulation matches exactly with the lattice. Beam-beam interaction at each the 
IP was represented by a kick in transverse phase space and the kick was calculated by using the 
particle-in-cell method as described in Ref. Since the beams in the HERA experiment were 
flat, a uniform mesh extending to ±200" in the configuration space with a grid constant of 0.2a 
was necessary in this case. Tracking of particle motion was conducted in four-dimensional trans- 
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verse phase space without synchrotron oscillations and momentum deviations. For lepton beams, 
the quantum excitation and synchrotron damping were treated as kicks in each turn during the 
tracking. The horizontal kick is ^S] 



where x and Px are the normalized horizontal coordinate and its conjugate momentum, Cx is the 
horizontal emittance, and wi and W2 are random numbers with a Gaussian distribution that is 
centered at zero and has unit standard deviation. The damping time in horizontal and vertical 
direction, Tx and Ty, have the unit of turns. For HERA, Tx = 436 and Ty = 600, respectively. The 
vertical kick has a similar formula. 

With the beam-beam simulation, the emittance growth of the e"*" beam and the specific lumi- 
nosity were calculated as functions of (,e,y for the HERA experiment and the results were plotted 
in Figs. 1 and 2 as solid lines. Both the emittance and the luminosity plot show a remarkable 
agreement between the experiment and the simulation. Figure 3 plots the evolution of the vertical 
emittance of the e"*" beam at different f3e,y In all these cases, after a quick emittance blowup, the 
beam emittance is re-stabilized and, consequently, an equilibrium (or quasi-equilibrium) state of 
the e"*" beam was reached. During the emittance blowup, the particles in the beam core escape 
quickly to the beam tails. Without the onset of the coherent beam-beam instability, the particles 
in the beam tails are stable for the beam-beam interaction. This re-stabilization of the beam emit- 
tance is therefore due to a depopulation of the beam core and is consistent with the experimental 
observation that the beam lifetime and operation conditions were good during the experiment even 
in the case of S,e,y = 0.272 jT6] . 

In the beam simulation, the coherent beam-beam tune shift was also calculated. Fig. 4 is 
the calculated power spectrum of the coherent oscillation of the e"*" beam at /3e,j/ = 4.0 m. Due 
to the quantum fluctuation, the e"*" beam has always a very small oscillation which is enough for 
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the calculation of the coherent frequency if the numerical noise is small. Since a large number of 
particles were used for each of the beams, in this simulation the numerical noise was very low and 
we were able to calculate the coherent frequency without applying off-center kicks on the beams. 
The spectrum was calculated by the fast Fourier transformation (FFT) of the beam centroid motion 
from the 5000th to 9000th turn. Since at the 5000th turn the beam has already reached its quasi- 
equilibrium, the transient state of the beam was thrown out during the tune calculation (see Fig. 
3). As shown in Fig. 4 the power spectrum peaks at 0.1605 and 0.2331 in the horizontal and 
vertical planes, respectively, which corresponds to Ux = 52.161 and Uy = 52.233 for the beam 
coherent tunes during collision. This simulation result agrees excellently with the experimental 
measurement of Vx = 52.160 and Vy = 52.233. The beam-beam simulation therefore confirmed the 
coherent beam-beam tune shift measured in the HERA beam experiment. 

To understand the very small coherent beam-beam tune shift in this case of very large beam- 
beam parameter, the dynamics of the beam particle distributions was studied during the beam 
simulation. As shown in Fig. 5, the distribution of the e+ beam deviates from a Gaussian distri- 
bution with a significant drop at the beam core and a growth of the beam tails. Compared with 
the distribution of the e"*" beam, a Gaussian beam that has the same emittance of the e"*" beam has 
more particles in the core. The real coherent beam-beam tune shift measured in the experiment and 
calculated with the beam simulation is therefore smaller than that calculated from the rigid-beam 
model with Gaussian beams. Moreover, the assumption of Gaussian beams in all the theoretical 
models of the coherent beam-beam oscillation neglects the effect of the (^dependence of the beam 
equilibrium distributions, where is the angle of the action-angle variable. In the HERA beam 
experiment, the equilibrium (quasi-equilibrium) distribution of the e+ beam in fact has a strong 
(?!)-dependence as shown in Fig. 6. For a comparison, the initial distribution used in the simulation 
was also plotted in the Figure and shows no ^-dependence as it should be. In order to compare 
the experimental/simulation result with calculations of theoretical models, the theoretical models 
need to be modified to include the (^dependence of non-Gaussian equilibrium distributions. 
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IV. VALIDITY OF THEORETICAL MODELS FOR BEAM COHERENT 

OSCILLATION 

Two theoretical models, the rigid-beam model and the linearized Vlasov equation with the 
single-mode approximation, were examined with the HERA experiment. For the rigid-beam model, 
the coherent tunes were calculated with two different methods: assuming a Gaussian distribution as 
the equilibrium distribution of the e"*" beam [Eqs. (AlO) and (A14) in Appendix A.l] or using the 
non-Gaussian quasi-equilibrium distribution of the e"*" beam obtained in the beam-beam simulation 
(Appendix A. 2). In the case of the Gaussian distribution, the beam sizes used in Eq. (A14) are 
of the experimental measurement or the beam-beam simulation. The small differences in the 
beam sizes measured in the experiment or calculated from the simulation (see Fig. 1) made little 
difference in Eq. (A14). For the linearized Vlasov equation, the horizontal coherent tunes were 
obtained by solving the initial- value problem of the linearized Vlasov equation with the single-mode 
approximation in the horizontal plane. The details of the calculations are in Appendix A and B. 

Tables II and III list the ratio of the coherent and incoherent beam-beam tune shifts of the e^" 
beam, A.V(.^x/ {'^ie,x) and Afe,j//(2^e,j/)) calculated in the beam simulation or calculated with the 
theoretical models for the cases of I3e,y = 1.0 and 4.0 m, where Afe^i and ^i'e,y are the horizontal 
and vertical coherent beam-beam tune shift of the beam. The experimental measurement 
at (5e,y = 4.0 m is also included in Table II. The significant discrepancy between the results of 
the models and the results of the experiment/simulation shows that the theoretical models are 
inconsistent with the experiment and simulation. Note that the rigid-beam model with the beam 
distribution from the simulation did a little better than that with Gaussian distribution. To further 
examine the failure of the theoretical models, a beam-beam simulation was conducted for the case 
of I3e,y = 1.0 m but with only one tenth of the p-bunch current used in the experiment, i.e. 
Ce,2/ = 0.0068. The result is listed in Table IV and shows that the rigid-beam model with either 
the Gaussian distribution or the distribution obtained from the simulation are in a good agreement 
with the beam simulation. In fact, the quasi-equilibrium distribution of the beam in this case 
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is very close to a Gaussian. This re-confirms many previous studies that the rigid-beam model is 
correct for unsymmetrical beam-beam interactions with relatively small beam-beam parameter. 

In the case of the linearized Vlasov equation, as shown in Tables II-IV, the calculation yielded 
wrong results no matter how small the beam-beam parameter is. In the calculation with the 
Vlasov equation, several approximations were employed. Among them, the linearization, the one- 
dimensional beam, and the single-mode approximation are the three major approximations that 
cannot be directly justified by the experimental observations (see Appendix B). The linearization of 
the Vlasov equation should not play the leading role in its failure for the unsymmetrical beam-beam 
interaction since the linearization should be valid for a relatively small beam-beam parameter. To 
verify the validity of the one-dimensional approximation, we did a series of beam-beam simulations 
for the case of Pe,y = 1-0 m, ^e,y = 0.068 or 0.0068, but with a different beam aspect ratio {cry^e/(^x,e) 
ranging from 0.063 to 0.120. In all these cases, the horizontal coherent beam-beam tune shift was 
found to be very similar to that of the HERA experiment. Note that in the HERA experiment 
and in all the cases in Tables II-IV, cry^e/o'x,e = 0.126. This indicates that the horizontal coherent 
beam-beam tune shift is not sensitive to the aspect ratio of the beam. Similar phenomenon has 
also been observed in the symmetrical case of beam-beam interactions 0. The large discrepancy in 
the coherent beam-beam tune shifts calculated with the linearized Vlasov equation as compared to 
that obtained from the experiment /simulation is apparently not due to the approximation of one- 
dimensional beam. With the single-mode approximation, on the other hand, the only oscillation 
mode of the beam distribution that was kept in the calculation is the m = 1 mode in Eq. (B6). 
As shown in Fig. 6, the equilibrium distribution of the beam [fi^Q{I,(f)) in Eq. (B4)] depends 
strongly on (j). In this case, the dominant mode of the equilibrium distribution is the m = 2 mode 
which has a ^-dependence of e'^'^'^. As the m = 2 mode could be a dominant intrinsic mode of 
the system, the single-mode approximation could be inconsistent with the dynamics of the beam. 
In order to study the effect of mode couplings in the linearized Vlasov equation, we have derived 
the eigenvalue equation of the linearized Vlasov equation that is similar to Eq. (B13) but includes 
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high-order modes. After truncating higher-order modes at m = rrip, M in Eq. (B13) becomes a 
(2mp -|- l)(2/p -|- 2) X (2mp + l)(2/p -|- 2) matrix where Ip is the number of grids on the mesh of the 
action space and Mi and M2 in Eq. (B14) are no longer diagonal matrices (see Appendix B). We 
have, however, failed in obtaining a set of orthogonal eigenvectors for the eigenvalue equation of the 
linearized Vlasov equation because M is a singular (ill-conditioned) matrix when the mode couplings 
are included. Similar problem has also been encountered when including the mode coupling in the 
case of symmetrical beam-beam interactions. 

V. COHERENT VERSUS INCOHERENT BEAM-BEAM TUNE SHIFT 

In the symmetrical case of beam-beam interactions, the ratio of the coherent and incoherent 
beam-beam tune shift is approximately a constant of 1.2 for round beams and 1.3 for flat beams 
[1]. Note that with symmetrical beam-beam interactions, the beam-size growths of the two beams 
are symmetrical and the beam distributions are usually close to a Gaussian when the beam-beam 
parameter is below the threshold of the onset of coherent beam-beam instability [21 El- The 
coherent beam-beam tune shift depends linearly on the beam-beam parameter in a fairly large 
range of beam-beam parameter On the contrary, in the HERA experiment the ratio of the 
coherent and incoherent beam-beam tune shift of the e"*" beam decreases monotonically with the 
increase of the beam-beam parameter as shown in Fig. 7. This different characteristic of the 
coherent beam-beam tune shift stems mainly from the mismatch in the equilibrium distributions 
of two unsymmetrical colliding beams. In the HERA experiment, the beam size of the e"*" beam 
at the IPs is slightly larger than that of the p beam initially without collision (see Table I). When 
the beam-beam parameter is small such as in the case of ^e,y = 0.0068, the beam-size growth is 
insignificant and the beams are close to the Gaussians. The ratio of the coherent and the incoherent 
beam-beam tune shift of the beam is simply determined by the mismatch in the initial beam 
sizes without collision as described by Eq. (A14) of the rigid-beam model. When the beam-beam 
parameter is large, on the other hand, the beam-size growth of the e"*" beam is significant (see Fig. 
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1) and dominates the beam-size mismatch. Moreover, the particle distribution of the e"*" beam 
deviates significantly from a Gaussian distribution (see Fig. 5). In this situation, the size of the 
e"*" beam is much larger than that of the p beam and a large number of e"*" particles escape to the 
beam tails that is away from the core of the p beam. The coherent beam-beam tune shift of the e+ 
beam thus becomes much smaller than that in the case of two matched Gaussian beams. Note that 
for two matched Gaussian beams, the ratio of the coherent and incoherent beam-beam tune shift is 
0.5 for unsymmetrical beam-beam interactions. As the beam-size growth and the resultant beam- 
size mismatch increase with the beam-beam parameter, the ratio of the coherent and incoherent 
beam-beam tune shift of the beam decreases with the beam-beam parameter. In general, the 
equilibrium states of the Vlasov equation in Eq. (B4) are different for two unsymmetrical colliding 
beams due to the asymmetry of beam-beam perturbations. The mismatch in the distributions as 
well as the beam sizes is, in principle, independent of the initial states of the beam distributions and 
the initial beam-size mismatch as long as the considered equilibrium state of the Vlasov equation 
in Eq. (B4) is well isolated and is the only one that is close to the initial beam distributions. The 
beam-size mismatch in the unsymmetrical case of beam-beam interactions is therefore intrinsic 
and unavoidable especially when the beam-beam parameter is larger. Moreover, this beam-size 
mismatch increases with the strength of beam-beam perturbations. The functional dependence of 
the coherent beam-beam tune shift on the beam-beam parameter in Fig. 7 is therefore a general 
characteristic of unsymmetrical beam-beam interactions. 

Figure 7 also plots l^Ve,yl {^Ce,y) as a function of ^e,y calculated based on the rigid-beam model. 
It shows that when the beam-beam parameter is large, the rigid-beam model results in a similar 
^e,j/-dependence of At'e,2y/(2^e,y) even though it overestimates the coherent beam-beam tune shift. 
Note that in the rigid-beam model, this ^e,j/-dependence is due to the beam-size mismatch as given 
in Eq. (A14). At a large ^e,y, the beam-size mismatch is dominated by the beam-size growth of 
the e"'" beam. This suggests that the result of the rigid-beam model in Eq. (A14) provides an 
approximation of the functional dependence of the coherent beam-beam tune shift to the intrinsic 
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beam-size mismatch in misymmetrical beam-beam interactions with large beam-beam parameters. 

VI. CONCLUSION 



The coherent beam-beam tune shift was studied in the case of unsymmetrical beam-beam 
interactions where the two beams have very different beam-beam parameters and betatron tunes. 
The results of a self-consistent beam-beam simulation, the rigid-beam model, and the linearized 
Vlasov equation were compared with the beam measurement in the HERA 2000 Beam Study. 
Remarkable agreement was found between the beam simulation and the HERA experiment in a 
wide range of and, especially, at very large beam-beam parameters of the lepton beam. The rigid- 
beam model was found to be only correct when the beam-beam parameter is small. The result of 
the linearized Vlasov equation with the single-mode approximation is inconsistent with the result 
of the beam experiment/simulation in either cases of large or small beam-beam parameter. The 
failure of the linearized Vlasov equation could be due to the single-mode approximation used in 
solving the linearized Vlasov equation. A study of the dynamics of the beam distribution showed 
that the high-order modes are important to the beam dynamics in this situation. An attempt 
to include high-order modes in the calculation has, however, not been successful because of the 
difficulty in finding a set of orthogonal eigenvectors for the linearized Vlasov equation. Recently, 
efforts have been made to include the angle-dependence of beam distributions in an expansion of 
the Vlasov equation (19) . More studies are needed for a relevant solution of the linearized Vlasov 
equation for the unsymmetrical beam-beam interaction. Currently, the numerical simulation is the 
only reliable approach for an prediction of the coherent beam-beam tune shift in this situation. 

One interesting phenomenon observed in this study is the very small coherent beam-beam tune 
shift in this unsymmetrical case of beam-beam interactions. It was found that the ratio of the 
coherent and incoherent beam-beam tune shift of the weak lepton beam in HERA decreases from 
0.3 to 0.02 as the incoherent beam-beam tune shift increases from 0.01 to 0.54. On the contrary, 
in the symmetrical case of beam-beam interactions, this ratio maintains approximately a constant 
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of 1.2 for a round beam or 1.3 for a flat beam in a large range of beam-beam parameter. The 
reason for this different characteristic of the coherent beam-beam tune shift is the intrinsic beam- 
size mismatch between two unsymmetrical colhding beams due to the difference in the equihbrium 
distributions of the two beams. This intrinsic mismatch in the beam distributions due to beam-beam 
interactions becomes more pronounced as the strength of beam-beam perturbations increases. The 
ratio of the coherent and incoherent beam-beam tune shift decreases, in general, with the increase 
of beam-beam parameters in the unsymmetrical case of beam-beam interactions. 

APPENDIX A. COHERENT BEAM-BEAM TUNE SHIFT FROM RIGID-BEAM 

MODEL 

The coherent beam-beam tune shift in the rigid-beam model has been studied before for the 
unsymmetrical beam-beam interaction. Previous studies are, however, limited to the special cases 
where either the two beams have the same or very different jl4j beam-beam parameters. More- 
over, all of them assumed that the beams are Gaussian. In the following, a general formula of the 
rigid-beam model is derived that can be apply to any case of unsymmetrical beam-beam interaction 
with either Gaussian or non-Gaussian beams. 

Let Pi (r, 9) be the distribution of beam i in the normalized configuration space where i = 1 
or 2, r = are the normalized coordinates of the transverse space, and 9 the azimuthal angle 

associated with the path length along the closed orbit. The beam centroid in the normalized space 
can be calculated by i?i = / fpi (r, 9) dr. Considering a linear lattice with one IP, the transverse 
motion of the beam centroid can be described by 



are the fractional parts of the betatron tunes of the lattice for beam i. The main part of the 





In Eq. (Al), rij is a 2 X 2 diagonal matrix with (r2j)ii = !/■ 



ix and {ni)22 



'ly, where {i^i^^, i^i,y) 
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beam-beam kick in Eq. (Al) is 

F= J J p, {ri,e) p2 (r2, e) G (/^J/.^i - 0^/^X2 , Pl^Jyi - (3l{^y2) dndr2 {A2) 

— oo — oo 

where G{x,y) = rjr^ is the Green Function of the beam-beam interaction and [fii^xi f^i,y) are 
the horizontal and vertical beta-function of beam i at the IP, respectively. The strength of the 
beam-beam kick for the horizontal component of Eq. (Al) is 

_ .1/2 

— M,x — '^i,xPix 

where i =1 or 2, j =1 or 2, but i 7^ j; Nj is the number of particles per bunch of beam j; ai the 
classical radius of the particle in beam i; and 7^ the Lorentz factor of beam i. With the definition 
of the beam-beam parameters in Eq. (1), this kick strength can be written as 

^i,x — -TJ2 (^j,x\(^j,x -r CFj^y) (^/lo; 
Pi,x 

where (^i,^) ^i,y) and ((Jj^^, Ui^y) are the beam-beam parameters and the rms beam sizes at the IP 
of beam i. For the vertical component of Eq. (Al), Aj = Xi^y can be easily obtained by exchanging 
X and y in Eq. (A3). 

In the rigid-beam model of the coherent beam-beam oscillation, the shapes of the particle 
distributions in phase space are assumed to not change with time during the beam oscillation while 
the centers of the distributions oscillate with the beams' coherent tunes. The distribution during 
the beam oscillation is thus assumed to be pi{r,9) = poi{f— Ri), where poi{r) is the equilibrium 
distribution when the beam is centered at the closed orbit. In general, F in Eq. (A2) is a function 
of moments of phase-space variables and the time-dependence of F is implicitly through all the 
moments. (To have a better picture of this, one may consider a moment expansion of beam particle 
distributions in phase space). With the rigid-beam approximation, F depends on the lowest-order 

moments Ri and R2 only, 

+00 +00 

F= J J poi (ri - Ri) P02 (r2 - R2) G ((i\{^xi - (3l{^X2 , (i\[yyi - f3l[yy2) drM {AA) 
—00 —00 
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Note that the rigid-beam model may fail in cases where the variation of the distributions is im- 
portant during the beam coherent oscillation. To find the oscillation frequencies of one can 
average the beam-beam kick in Eq. (Al) over one turn (27r in longitudinal direction) and expand 
F into a Taylor series of Ri. Keeping only the linear terms of Ri, Eq. (Al) becomes a coupled 
four-dimensional harmonic oscillator, 



where i = 1 or 2, and Aj are 2x2 matrices with 

dF 

Ai = 



ni-Ri = i-iy+^Xi (Ai ■Ri+A2- R2) {A5) 



If both the beams are mirror symmetric with respect to the horizontal and vertical plane, Aj 
are diagonal matrices and the horizontal and vertical coherent oscillation are decoupled. The two 
eigenfrequencies for the coherent oscillation in the horizontal plane can then be solved as 



where 



{A8) 



Suji = Xi^x (Ai)ii / {2vi^x) 
for i = 1 or 2. If 5u)i « Ui^x and in Eq. (A7) 



{uJl - ul)"^ » lQui^x^2,x50Ji5uJ2 (^9) 

then the coherent frequencies of the two beam are 



(^10) 



V- =L02 = l'2,x - 5UJ2 

where and are the horizontal coherent tunes of beam 1 and beam 2, respectively. Note that 
in the unsymmetrical case of beam-beam interactions, the two eigenfrequencies do not correspond 
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to the so called 0- (or a-) and 7r-modes of symmetrical beam-beam interactions. The condition in 

Eq. (A9) can be further simplified as 

Wi,x - ^2,x\ » VSuJiSuJ^ ■ (^11) 

Therefore, if the difference of the lattice tunes is much larger than the geometric average of the 
coherent beam-beam tune shifts of the two beams, the coherent beam-beam tune shifts can simply 
be calculated with Eq. (AlO). Note that in the HERA experiment this condition was fulfilled. For 

the case of strong- weak beam-beam interactions such as ^i_x >> (,2,x, one can expand iy± in terms 
of ^2,x/^i,x- Keeping only the dominant term in the coherent beam-beam tune shifts yields 

Z^— = 1^2 T 2 2 TT r 0U2 

where z/_|_ and are the coherent frequencies of the weak (beam 1) and strong (beam 2) beam, 
respectively. In the first equation of Eq. (A12), since the zeroth-order term {5ui\) of \2,x/^i,x exists 
and dominates the coherent beam-beam tune shift of beam 1, the first- or higher-order terms were 
neglected. In the second equation of Eq. (A12), on the other hand, the zeroth-order term is zero and 
the first-order term was thus kept. Note that if the lattice tunes of the two beams are very different, 
Eq. (A12) is equivalent to Eq. (AlO). If the denominator in Eq. (A12), z^^^, — ff ,,. — 2vi^xSu)\, is 
small, an analysis of the higher order terms shows that the expansion in terms of X2,x/^i,x is no 
longer accurate and the coherent tunes have to be calculated by using Eq. (A7) . The two coherent 
frequencies in the vertical plane can be easily obtained by exchanging x and y and changing (Aj)j^^ 
with (Aj)22 in Eqs. (A7)-(A12). As shown in our study (see Section IV), this approach of the 
rigid-beam model is quite good in the case of unsymmetrical beam-beam interactions with a small 
beam-beam parameter. 
A.l. Gaussian Beams 

In the case that poi are Gaussian distributions, matrix Aj in Eq. (A6) can be calculated 
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analytically with Eq. (A4) as 



(A.) 



11 ~ 




{A13) 



where Ti^ = 



'y ~ 




The matrix element of (Aj)22 can be simply 



obtained by exchanging x and y in Eq. (A13). Substituting Eqs. (A2) and (A13) into Eq. (A8) 
yields 



where i = 1 or 2; j = 1 or 2; but j 7^ i. When i = 1, Eq. (A14) gives the coherent beam-beam 
tune shifts of the weak beam [beam 1 in Eq. (A12)] obtained previously by Hoffstaetter for the 
strong-weak case of beam-beam interactions |14) . If {vi^x-,i^i,y) = i'^2,x , '^2,y) and {cri^x,o'i,y) = 
{cr2,x, (^2,y), Eq. (A7) is reduced to the formula obtained by Hirata If {S,i,xi(,i,y) = {S,2,x,(,2,y) 
and {cri^x,cri,y) = {(^2,x, '^2,y), Eq. (A7) is reduced to the formula obtained by Hofmann |H]. 
A. 2. Non-Gaussian Beams 

For non-Gaussian beams, especially the distributions obtained from beam-beam simulations 
such as that in Figs. 5 and 6, matrices Aj in Eq. (A6) cannot be obtained analytically but can be 
calculated numerically by using Eq. (A4). The coherent frequencies can then be calculated with 
Eq. (A7) or directly from Eq. (A5) if the horizontal and vertical motion are coupled. 



APPENDIX B. COHERENT BEAM-BEAM TUNE SHIFT FROM LINEARIZED 



The use of the linearized Vlasov equation has been very successful for the coherent beam- 
beam tune shift in the case that two beams have the same or very close lattice tunes 111. In 
order to find the coherent beam-beam tune shift, one needs to identify the coherent frequencies 
from the eigenfrequencies of the linearized Vlasov equation. The linearized Vlasov equation, in 
principle, has infinite numbers of eigenfrequencies associated with infinite numbers of oscillation 
modes. For real beams, the number of the eigenfrequencies of the beam oscillation in transverse 
space is twice of the number of particles in a bunch. When two beams have the same lattice tunes. 



.l.X 




(AU) 



VLASOV EQUATION 
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the coherent frequencies can be easily identified since the eigenfrequencies that correspond to the 
coherent frequencies are separated from the rest of the eigenfrequencies that form a continuous 
band (many close eigenfrequencies lines) 1^. Although, it is not very clear mathematically why this 
separation occurs. The situation becomes more complicated when two beams have very different 
lattice tunes. In this case, all the eigenfrequencies are in one or two continuous bands and the 
coherent frequencies cannot be identified by only solving the eigenfrequencies. In order to find the 
coherent frequencies in the HERA beam experiment, we will instead solve the initial- value problem 
of the linearized Vlasov equation for the coherent beam oscillation. 

Consider only the horizontal motion (very flat beam) in a linear lattice with one IP. In terms 
of action-angle variable, the Hamiltonian for the betatron motion of beam i (i = 1 or 2) can be 
written as 



Hi{i,(i>,e) = Hi,o{i) + u,ii,<i),e) 



27r^5(6' - 2irn) 



(Bl) 



where Hi^ = Vi^xl is the Hamiltonian associated with the betatron motion in the linear lattice and 
Ui is the potential energy for the beam-beam interaction that can be written, for one-dimensional 
beams, as 

2 27r oo 



Ui{I,4>,e) = Uiifj] = -2^^^-^ J J fj{I'A',d)ln (y'2A,,/sin0 - ^2/?,- ,1' sin 0' j {B2) 



where i = 1 or 2, and j = 1 or 2, but i ^ j- The action-angle variables are related to the normalized 
variables by x = \/2Jsm(f) and p = V2Icos((>. fi{I,<p,9) is the particle distribution of beam i in 
phase space and satisfies the Vlasov equation. For convenience, we also define a functional Ui[fj] in 
Eq. (B2) for the potential integral. In Eq. (B2), ln(x — x') is the Green function for the potential 
of beam-beam interaction in one-dimensional space. If only the coherent beam-beam tune shifts 
are interested, one can get rid of the periodic 5-function in the Hamiltonian in Eq. (Bl) by average 
the beam-beam force over one turn. The Vlasov equation for /j can then be written as 
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where { } is the Poisson bracket. Assume that the beams have reached equihbrium distributions 
fifi that satisfy 

= {U^,o, ko} (54) 

where Uifi{I,(j)) = Ui[fjfi]. Consider that beam i experiences a smah perturbation from its equi- 
hbrium distribution 'ipi{I,(j),9) = fi{I,(j),9) — fifi{I,(j)). The hnearized equation for ^i{I,(j),6) can 
be obtained by subtracting Eq.(B4) from Eq. (B3) and neglecting the term {Uili/jj], ipi} which is 
higher-order in -0^ as 

where V;(/,(/.,0) = C/i[^,]. 

To solve Eq. (B5), one can convert it into a system of infinite numbers of coupled ordinary 
differential equations of modes by using Fourier transformation 

y oo 

-oo "^=-00 

where v is the oscillation frequency of the beams and m numbers modes. The m = 1 mode 
corresponds to the coherent dipole oscillation. To further simplify the problem, one may use the 
single-mode approximation in which only the mode with m = 1 is kept in the linearized Vlasov 
equation [3]. It turns out that the use of single- mode approximation is not only a convenience but 
also a necessity. Without the single-mode approximation, no effective method is available for the 
general solution of Eq. (B5), except for simplified models as in Ref. jTHI. Substituting Eq. (B6) 
into Eq. (B5), multiplying e"*"^ and integrating (f) over 27r on the both sides of Eq. (B5), and only 
keeping the m = 1 mode yields 





where = 
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and 



21 cos (f> 



dfifl 



— sm 




7r^2/5j,^/ J y/Pi^r,! sin - y/(3j,xl' sin 0' 



a/ 



If the equilibrium distributions are independent of (j) such as for Gaussian beams, the imaginary term 
of Gi{I,I') is zero. Otherwise, this imaginary term contributes a damping to the hnearized Vlasov 
equation when tjji is stable or an excitation when ipi is unstable. If the equilibrium distributions 
fifi are Gaussian, with a similar algebraic treatment in Ref. ;4, the integrals in Eqs. (B8) and 
(B9) can be calculated analytically as 



where i = 1 or 2; j = 1 or 2; but i / j. Zi = Pi^xl/crf ,^, z'j = Pj^^I' / ^'j^, and rij = (j'j^^/af^^. 

B.l. Eigenfrequencies and Eigenvectors of Linearized Vlasov Equation 

To further proceed with Eq. (B7), one may discretize the action space (I) into a mesh and 
solve the equation on the grids [311]. Let / = /A/ where A/ is the grid size; / = 0, 1, 2, ... Ip] 
and /pA/ is the size of the mesh. Since the distributions decay to zero quickly as / increases, a 
mesh that covers several (Ji^x is good enough for a calculation of the coherent frequency. In order 
to have an accurate frequency for the lattice tune in the eigenfrequencies of Eq. (B7), however, the 
mesh has to be large enough so that the incoherent beam-beam tune shift at / = lpA.1 is negligible. 
In this study, we therefore used lpIS.1 = 160 ei^x and A/ = 0.05 ej^^i where ei^x is the normalized 
emittance of beam i. Let '4>i{lAI,i') = Eq. (B7) can then be converted into a system of 

linear algebraic equations on the mesh, 




(SIO) 




i^^ii = i^i,xAi + Qi A/) + A/ ^ Gi (/A/, kAI) il^jk ■ 



(512) 



k=0 



that leads to an eigenvalue problem 



MV = uV 



(513) 
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where 



V = (Vio, ^11, • • •, V^i/p, ■020, -021 



2ii 



M is a 2(/p + 1) X 2(/p + 1) matrix 



M 



^ Ml Oi ^ 



V 



(514) 



O2 M2 

where Mi, M2, Oi, and O2 are (/p+1) x (/p+1) matrices. Because of the single-mode approximation, 
Ml and M2 are diagonal matrices with the diagonal elements 



(M. 



and the elements of Oj are 



(0.)h = G, {{k - 1)A/, (/ - 1)A/) 



{B15) 



(516) 



where k = 1, {Ip + 1); / = 1, {Ip + 1); and z = 1 or 2. If the equilibrium distributions are 
Gaussian, all these matrix elements in Eqs. (B15) and (B16) can be calculated analytically by 
using Eqs. (BIO) and (Bll). In the case of the HERA experiment, the equilibrium distribution of 
the p beam is still very close to a Gaussian but the e+ beam is no longer a Gaussian beam (see 
Figs. 5 and 6). Let beam 1 and 2 be the and p beam, respectively. Mi and O2 can then be 
obtained analytically. The matrix elements of M2 and Oi, on the other hand, have to be calculated 
numerically by using Eqs. (B8) and (B9) with the quasi-equilibrium distribution of the beam 
obtained from the beam simulation. 

With the eigenvalue equation in Eq. (B13), the eigenfrequencies and a set of orthogonal eigen- 
vectors for the linearized Vlasov equation can be found numerically. If the two beams have the 
same lattice tune, the eigenfrequencies of Eq. (B13) are identical to that obtained in Ref. 1^. 
Fig. 8a is an example of eigenfrequencies of Eq. (B13) for the case of i^i^x = '^2,x and ^i^a; = 6.2,x- 
It shows that in the symmetrical case of beam-beam interactions the coherent frequency (the first 
frequency line from the left of Fig. 8a) is separated from the rest of the eigenfrequencies that form a 
continuous band. The width of the band equals the incoherent beam-beam tune shift. The coherent 
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frequency can therefore be easily identified in this case. Note that the coherent beam-beam tune 
shift calculated from this coherent frequency is the same as that in Ref. [1]. The situation is more 
complicated when the two beams have very different lattice tunes. Fig. 8b plots the eigenfrequen- 
cies for the case of the HERA experiment. In this case, the eigenfrequencies are divided into two 
groups, one for each beam. For the e"*" beam, the eigenfrequencies form a continuous band that 
starts at the lattice tune of the e"*" beam and has a width of the incoherent beam-beam tune shift 
of the e"*" beam. Because there were two IPs in the HERA experiment, the incoherent beam-beam 
tune shift in Fig. 8b is 2^e,x = 0.082 for the e~*" beam. The characteristics of the eigenfrequencies for 
the p beam, in principle, is similar to that of the beam. Since the beam-beam parameter of the 
proton beam is very small {^2,x = ^p,x ~ 10~^), all the eigenfrequencies for the p beam degenerate 
into a single line (the first line from the left of Fig. 8b) that corresponds to the lattice tune of 
the p beam. In the case of very unsymmetrical beam-beam interactions, therefore, the coherent 
frequencies cannot be simply identified from the eigenfrequencies of the linearized Vlasov equation. 

It should be noted that because of a very small beam-beam parameter of the p beam, during 
the HERA experiment no coherent beam-beam tune shift was observed on the proton beam. The 
matrix elements of O2 are very small as compared with the diagonal elements of Mi and M2. 
O2 can therefore be approximated as a zero matrix and the eigenfrequencies for the e"*" beam can 
be easily obtained from MiVi = vVi, where Vi = (^ipio, ipu, ■ ■ •, V^i/p) is the sub vector space 
associated with the e"*" beam. Since Mi is diagonalized, solving the eigenfrequencies and a set of 
orthogonal eigenvectors of Mi is trivial. The eigenfrequencies and eigenvectors obtained from Mi 
were found to be the same as that of Eq. (B13) in the sub vector space associated with the e"*" 
beam in this case. 

B.2. Initial- Value Problem for Coherent Frequencies 

Let (i/i, • • •, z.,^+1, ■ ■ ; ^^2/^+2) and {v^, ■ ■ •, yC^+D, V('^+'\ ■ ■ •, y(2'p+2)) be the 

eigenfrequencies and eigenvectors of the discretized and linearized Vlasov equation where MF^"^ = 
f„y'("'). In the HERA experiment, the lattice tunes of the two beams are very different and, 
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therefore, the two eigenfrequency bands of Eq. (B13) are well separated (see Fig. 8b). In this case, 
{f„} are the eigenfrequencies for the e+ beam when n = l,...,lp + 1 and the eigenfrequencies for 
the p beam when n = Ip + 2, ...,2lp + 2. In the discretized action space, the perturbation of the 
beam distribution 9) and V'2(-f, ^, 0) can be represented as a vector. 

With the single-mode approximation, the general solution of <p, 9) and 'ip2{I^ <P, 9) can then be 
obtained from a superposition of the eigenvectors of the linearized Vlasov equation, 

2/p+2 

^{e) = CkV^^^e'^^-"'^'^^ {BIT) 

k=l 

where {Cfe} are constants and can be determined with an initial condition, (f), 0) and tp2^, 4>-> 0)- 
Since |Cfep is the oscillation amplitude of the beam distributions with the frequency of v^, the 
diagram of |Cfcp v.s. Vk corresponds to the frequency spectrum of the coherent oscillation. The 
two peaks in the \Ck\^-Vk diagram, therefore, provides the coherent tunes when ^, 0) — > and 

Consider a small kick that kicks beam 1 away from its equilibrium distribution fi,o{I, (p), where 
fifi{I,(j)) is known numerically from the beam-beam simulation. The initial perturbation of the 
beam distribution is 

^^1(1,0,0) = fi,o{x + xo,p^) - fi,o{x,p^) = Y,gm{I)e'""^ {B18) 

m 

and iIj2{I,<P,0) = where xq is the initial kick. With the single-mode approximation, , 0) ~ 
gi{I)e''^ and if{0) = g e*"^ where g = ( gi{0), gi{AI) gi{lpAI), ,...,0 )^. Note that the 
second half of the vector are all zero because beam 2 is not kicked. On the other hand, from Eq. 
(B17) 

(2lp+2 \ 
E CkVkje"^={'^c)e'^ (519) 

where V is a (2Zp-|-2) x (2Zp-|-2) matrix of which the ith column is Vi and C = (Ci,C2, C2/P+2) ■ 
The coefficients {C^} can then be calculated from C = \~^g. It should be noted that the initial 
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kick on the beam distributions in Eq. (B18) can be in any direction in phase space since the coherent 
frequency is the frequency of an infinitesimal oscillation. For near-integrable systems considered 
in this study, the phase-space region in the vicinity of the origin is integrable and only consists of 
invariant circles (tori). It is therefore isotropic. The coherent frequencies calculated were indeed 
found to be independent of the direction of the initial kick. 

Figure 8 plots the calculated |Cfep-ffe diagrams for the symmetrical case of beam-beam inter- 
actions where ui^x = i'2,x and ^1,3; = ^2,1 (Fig. 9a) and for the HERA experiment (Fig. 9b). In 
Fig. 9a, the peak with an arrow is the calculated coherent frequency that is the same as that in 
Fig. 8a. In Fig. 9b, the main peak indicates the calculated coherent frequency of the beam 
in the HERA experiment. The small peak in the lower right corner is the coherent frequency of 
the p beam. Figure 10 plots the calculated coherent beam-beam tune shift of the e"*" beam as 
a function of initial kick xq for the case of the HERA experiment. It shows that the calculated 
coherent beam-beam tune shift increases with the decrease of a;o and converges as xq approaches 
zero. This amplitude dependence of the coherent frequency is consistent with the beam simulation. 
Since the coherent frequency is the frequency of an infinitesimal oscillation, the convergence of the 
calculated coherent frequency at xq ^ provides the wanted coherent frequency. As shown in our 
study (see Section IV), the linearized Vlasov equation with single-mode approximation is a valid 
approach for the coherent beam oscillation with symmetrical beam-beam interactions but not with 
unsymmetrical beam-beam interactions. 
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Table I. Some beam parameters used in HERA 2000 beam experiment. The horizontal size of the 
e"^ beam at the IPs is (Te^x = 283 jjcm.. The horizontal and vertical size of the p beam at the IPs are 
(7p_a;=164 /xm and cre,2/=39.9 /xm, respectively. 



(^e,y (m) 


Ie+ (mA) 






^p,x/^p,y (10 ^) 


1.0 


19 


35.8 


0.041/0.068 


2.54/1.40 


1.5 


18 


43.8 


0.041/0.102 


2.35/1.06 


2.0 


17 


50.6 


0.041/0.136 


2.18/0.85 


3.0 


3.5 


62.0 


0.041/0.204 


0.43/0.14 


4.0 


2.6 


72.0 


0.041/0.272 


0.31/0.09 



Table II. The coherent beam-beam tune shifts of the e"*" beam at /3e,j/ = 4.0 m {^g^y = 0.272). 
"Experiment" and "Simulation" are the coherent tunes measured in the experiment and calculated 
in the beam simulation, respectively. "Rigid-Real" and "Rigid-Gaussian" are the coherent tunes 
calculated by using the rigid-beam model with a Gaussian distribution and with the distribution 
from the simulation, respectively. "Vlasov Eq." is the coherent tunes calculated with the linearized 
Vlasov equation. 



Pe,y = 4.0 m 


^e,x 


Az^e,x/(2^e,a:) 






Experiment 


0.1600 


0.110 


0.2330 


0.024 


Simulation 


0.1605 


0.104 


0.2331 


0.024 


Rigid-Real 


0.1555 


0.164 


0.2194 


0.049 


Rigid-Gauss 


0.1531 


0.194 


0.2040 


0.077 


Vlasov Eq. 


0.1123 


0.69 
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Table III. The same as Table II, but for /?e,y = 1.0 m {^e,y = 0.068). 



= 1.0 m 








Az^o ,,/ f2fe 


Simulation 


0.1600 


0.110 


0.2172 


0.212 


Rigid-Real 


0.1517 


0.212 


0.2074 


0.284 


Rigid-Gauss 


0.1475 


0.263 


0.1996 


0.341 


Vlasov Eq. 


0.121 


0.58 







Table IV. The same as Table III, but with only one tenth of the p-bunch current used in the 
experiment {^e,y = 0.0068). 



(3e,y = 1.0 m 




Az/e,x/(2Ce,a;) 


^e,y 


Aiye,y/{2ie,y) 


Simulation 


0.1672 


0.224 


0.2414 


0.337 


Rigid-Real 


0.1668 


0.263 


0.2406 


0.396 


Rigid-Gauss 


0.1668 


0.268 


0.2406 


0.399 


Vlasov Eq. 


0.165 


0.53 
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FIGURE CAPTIONS 

Figure 1. Emittance of the beam as a function of ^e,y is the emittance without colhsion. 
Discrete points are from the experiment and continuous curves from the beam simulation. Circles 
and curve a is the vertical emittance. Crosses and curve b are the horizontal emittance. The two 
experimental data points at each ^e,y where the measurement was performed correspond to the 
measurements at HI and ZEUS. 

Figure 2. The specific luminosity as a function of ^g.y Circles are from the experiment and 
continuous curves from the beam simulation. The two experimental data points at each ^e,j/ where 
the measurement was performed correspond to the measurements at HI and ZEUS. Crosses are the 
luminosity calculated with the measured emittance in Fig. 1 assuming Gaussian beam distributions. 

Figure 3. Evolution of the vertical emittance of the e+ beam calculated with the beam simulation 
for the cases of HERA beam experiment at (a) fi^^y = 1.0 m; (b) fi^^y = 1.5 m; (c) jd^^y = 2.0 m; (d) 
I3e,y = 3.0 m; and (e) (3e,y = 4.0 m. 

Figure 4. Power spectrum of the centroid motion of the beam in (a) horizontal and (b) 
vertical direction for the case of HERA beam experiment at ^e,y = 0.272 (/3e,j/ = 4.0 m). The beam 
centroid motion was calculated during the beam simulation. Note that in this case the incoherent 
beam-beam tune shifts for the two IPs are 0.082 and 0.544 in the horizontal and vertical plane, 
respectively. 

Figure 5. The projection of the particle distribution of the e+ beam in the vertical direction 
obtained by the beam simulation for the case of HERA beam experiment at ^e,y = 0.272 {I3e,y = 4.0 
m). (a) The initial Gaussian distribution, (b) the distribution at the 5000th turn, and (c) the 
Gaussian distribution of which the standard deviation is the same as that of the beam distribution 
in (b). 

Figure 6. The angle dependence of the particle distribution of the e+ beam obtained by the 
beam simulation for the case of HERA beam experiment at ^e,y = 0.272 (/3e,y = 4.0 m). (a) The 
initial Gaussian distribution, (b) the (jiaj-dependence and (c) the (;^j^-dependence of the distribution 
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at the 5000th turn. 

Figure 7. The ratio of the coherent and incoherent beam-beam tune shift as a function of the 
beam-beam parameter of the e"*" beam in the vertical plane calculated from (a) the beam-beam 
simulation and (b) the rigid-beam model with the Gaussian distribution. 

Figure 8. Eigenfrequencies of Eq. (B13) for the case of (a) i^i^x = t^2,x, ^i,x = ^2,x and (b) the 
HERA beam experiment. In (a), the dashed line marks the lattice tune that corresponds to the 
0-mode. The single isolated line on the left is the coherent frequency that corresponds to a ratio 
of the coherent and incoherent beam-beam tune shift of 1.35. In (b), the dashed lines mark the 
lattice tunes of the e"*" beam (the left line) and p beam (the right line), respectively. The single 
isolated line on the right is the degenerated eigenfrequencies for the p beam and the band on the 
left is for the e+ beam. The vertical axis has no physical meaning. 

Figure 9. |Cfep as a function of Uk for the cases of Fig. 8. The arrows indicate the coherent 
frequencies. 

Figure 10. Calculated coherent beam-beam tune shift by using the linearized Vlasov equation 
as a function of initial kick xq on the distribution of the e+ beam [see Eq. (B15)] for the case of 
the HERA experiment at /?e,y = 1.0 m. 



32 



Spectrum 



o 




p(y)/Po(o) 




27rp(0) 




